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ABSTRACT 

Stellar migration is an important dynamical process in Galactic disk. Here we model the radial stellar migration in the Galactic 
disk with an analytical method, then add it to detailed Galactic chemical evolution model to study the influence of radial stellar 
migration on the chemical evolution of the Milky Way, especially for the abundance gradients. We found that the radial stellar 
migration in the Galactic disk can make the profile of the G-dwarf metallicity distribution of the solar neighborhood taller and 
narrower, thus it becomes another solution to the "G-dwarf problem". It can also scatter the age-metallicity relation. However, 
after the migration, the abundance distributions along the Galactic radius don't change much, namely the abundance gradients 
would not be flattened by the radial stellar migration, which is different from the predictions of many theoretical works. But it 
can flatten the radial gradients of the mean chemical abundance of stars, and older stars possess flatter abundance gradients than 
younger stars. The most significant effect of the radial stellar migration on the chemical abundance is that at a position it scatters 
the abundance of stars there from a relatively concentrated value to a range. 

Subject headings: Galaxy: abundances - Galaxy: evolution -Galaxy: disk -Galaxy: kinematics and dynamics 



1. EVITRODUCTION 

The evolution of the Milky Way is clearly manifested by its 
structure formation and the enrichment of chemical elements 
heavier than those stemming from the primordial nucleosyn- 
thesis, and the latter resulting from the consumption of gas 
to star formation. Gas collapses to form stars, and various 
chemical elements are created by nucleosynthesis processes 
in stars. Through several ways, such as mass loss and super- 
nova explosion (during which heavier elements can also be 
produced), the newly synthesized nuclear products go back 
into the interstellar medium (ISM). Then the gas containing 
newly produced material forms stars of next generation, and 
the processes above will repeat again and again. The cycle 
between ISM and stars continuously enrich the Milky Way 
chemically on both the variety and the amount of chemical el- 
ements. Thus the dynamical behaviors of gas and stars would 
significantly influence the Galactic chemical evolution and 
contribute to form the present chemical features of the Milky 
Way. 

To understand how the Galaxy formed and have evolved 
on both structure and chemical composition, models of the 
chemical evolution of the Milky Way play an important 
role in the research. From the earliest closed-box Simple 
Model to the open models with gas in fall from outside of the 
Galaxy (Matteucci & Greggio' 1986", 'Matteucci & F rancoisI 
[1989; Chiappini et al. 1997; Portinari et al. 1998), and to 
the models including radial gas flow (Lacey & Fall Il985t 
IFortinari & Chios i 2000; Spitoni & Matteucci 2011, the first 
one doesn't relax IRA), the chemical evolution model of the 
Milky Way has been developing all the time. Up to now the 
effects of gas part on the Galactic chemical evolution has been 
studied in detail. 

However, the topics of how stars evolve dynamically in 
galactic disks and its effects on galactic chemical evolution 
have received less attention in the aforementioned works. One 
main reason resides that it has been believed for decades that 
stars would remain at roughly their birth radii throughout their 
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lives. Thus, according to the mechanism of chemical en- 
richment, there would be a perfect relation between the age 
and metallicity of stars (the so-called Age-Metallicity Rela- 
tion, AMR) at a certain position. For example, in the so- 
lar neighborhood the metallicity of stars decreases with the 
increasing of their age. This relation has been adopted by 
man y chemical evolution models in the last decades of years 
(e.gJTins\e^Wf3, 'Matteucci & Franco isl[T989tlCarigil l 19961; 
Chiappini et al. 1997; Boissier & Prantzosl ll999 ). As a con- 
sequence, the influences of the dynamical behaviors of stars 
on the Galactic chemical evolution have not been paid enough 
attention. 

However, as the enrichment of observational data and the 
development of theoretical studies in recent years, the classic 
picture has been replaced by new knowledge. Observation- 
ally, the tight age-metallicity relation has not been found in 
recent data. On the contrary the distribution of the age and 
metallicity of stars in the solar neighborho o d is flat together 
with large scatter (lEdvardsson et al.l [19931 iNordstrom et alj 
120041). and the scatter would increase with the increasing age 
jHavwoodl [20081; ICasagrande et all 1201 Ih . It suggests that 
stars born at other regions with diff'erent age and metallicities 
could appear in the solar neighborhood through migration. 
Besides, the metallicity distribution function (MDF) in the so- 
lar neighborhood, which acts as a constraint of the ch emical 
evolution model, obtained by ICasagrande et aTl (1201 lb shows 
that old stars have a considerably broader distribution than 
young stars while the peak always remains around the solar 
value, suggesting that the wings of MDF could mostly com- 
prise stars bom at various Galactocentric radii and migrated to 
the current position over d ifferent timescales. Theoretically, 
ISchonrich & BinnevI (l2009i) has distinguished two drivers of 
radial migration: by scattering at an orbital resonance or by 
non-resonant scattering by a molecular cloud, the angular mo- 
mentum of a star can be changed so that the star's guiding- 
center radius changes and the star's entire orbit moves inwards 
or outwards depending on whether its angular momentum is 
lost or gained, or the angular momentum of a star can be main- 
tained while the star's epicycle amplitude increases so that 
the star contributes to the density over a wider range of radii. 
They call them churning and blurring respectively. Besides, 
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iMinchev & FamaevI ( 120101) identified another radial migration 
mechanism in barred galaxies, namely the resonance overlap 
of the bar and spiral structure induces a nonlinear response 
leading to a strong redistribution of angular momentum in the 
disk. As a consequence, besides gas flows, the radial stellar 
migration is another important component that should be con- 
sidered in the Galactic chemical evolution model. 

When the radial stellar migration is included in chem- 
ical evolution models, the classical tight correlation be- 
tween age and metallicity is relaxed, and this relation be- 
comes eff'ectively an additional constraint independent from 
the metallicity distribution. Abundance gradients are also 
very important constraints for t he Galactic chemic a l evo- 
luti on model. As suggested by ISellwood & BinnevI ( 12002b 
and lSchonrich & BinnevI (l2009h . since the migration of stars 
makes it possible to transport materials to anyplace along the 
disk no matter where they originally located, the mixture of 
chemical elements over large range of radius would reduce 
the differences of abundances between different regions along 
the disk and flatten the abundance gradients. 

This paper is to study the influence of radial stellar migra- 
tion in the disk on the chemical evolution of the Milky Way 
and the degree to which it could affect the formation of some 
Galactic typical chemical features. We organize it as follows: 
in Section 2 we describe the reference models used in this 
work; the description of the radial stellar migration is pre- 
sented in Section 3; Section 4 and 5 show the stellar ejecta 
and the observational data adopted here, respectively; in Sec- 
tion 6 we report and discuss our results; and finally we draw 
our main conclusions in Section 7. 



2. GALACTIC CHEMICAL EVOLUTION MODEL 

The model adopted here is based on the model of 
iPortinari & Chio s? (7000), which is an open model where the 
disk forms gradually by accretion of protogalactic gas and 
radial gas flow exists. The disk is divided into N concen- 
tric rings or shells, each of which is of 1 kpc wide; and in 
each ring the gaseous component and its chemical abundances 
evolve due to: 1. Star formation which locks up gas into 
stars; 2. Infall of primordial protogalactic gas; 3. The ejecta 
of stars which were bom in the ring or migrated from other 
rings, bringing the enriched material back to the interstellar 
medium; 4. Gas exchange with the neighboring rings as a 
result of radial gas flows. 

The basic equation describing the variation of gaseous com- 
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where cr, (r, t) is the surface gas density of element i at radius r 
at time t, X,(r, t) - cTiir, t)la-gas{r, is the abundance by mass 
of element i, i/'(r, f) is the star formation rate at radius r at time 
t, 0(ot) is the initial mass function, /?„„(?) is the mass frac- 
tion of a star with mass m ejected into the ISM in the form 
of element i, and EsNebi is the mass contribution of a type la 
supernova explosion for element i. The first term represents 
the reduction caused by star formation. The following four 
integral terms describe the contributions of stars in different 
mass ranges. The summations represent the effects of stel- 
lar radial migration from all the considered radii, and fsM is 
the distribution function of radius describing the stellar migra- 
tion which will be shown in detail in Section 3. The last two 
terms represent the contributions of gas accretion and radial 
gas flow, respectively. 

2.1. Star formation rate 

The star formation rate (SFR ) [jj{r, t) adopted here is of the 
form in iTalbot & ArnetB ([l975h . 



cr{r, t) o-g(r, t) 
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where v is the efficiency of the star formation rate, cr^(r, t) 
is the surface mass density of gas, (x(r, f) is the total surface 
mass density, t) is the total surface mass density at a par- 
ticular Galactocentric distance r which is taken here to be 
T^o = 8.5 kpc. 

2.2. Initial mass function 

The initial mass function (IMF) we adopt is assumed to be 
constant in space and time, and the prescription for it are taken 
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from [Scaiol([T986h : 



<p(m)dtn — cm '^dm 



(3) 



where = 1.35 for m < 2 Mq and fi - l.l for m > 2 Mq. 
Since the majority of the chemical enrichment is due to stars 
of m > 1 Mq, it is meaningful to fix the fraction f of the total 
stellar mass distributed in stars above 1 Mq. Thus the lowest 
stellar mass considered here, M/ , is also fixed according to 
the normalization condition. 



(p(m)dm — 1 
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We adopt the upper mass limit, M,„ to be 120 Mq. 

2.3. Gas accretion 

An infall rate of primordial protogalactic gas decreases ex- 
ponentially in time with a timescale T(r): 



o-infirj) = A(r) e 



(4) 



where r(r) is chosen for two conditions: be constant at all 
radius and be variable with Galactic radius under an inside-out 
disk formation scenario. A(r) is obtained by reproducing the 
current total surface mass density distribution, and the present 
surface mass density of the disk is exponential with the scale 
length Rd -3.5 kpc and equals 54 Mq pc^^ at the solar radius, 
and 1 7 Mq pc^^ was adopted for the halo (IChiappini et al.l 
l200l . 

2.4. Stellar lifetime 

To drop the "instantaneous recycling approximation 
(IRA)", the finite lifetimes of stars of different masses must 
be considered. We adopt th e metallicity-dependent stellar 
lifetimes from iPortinari et all (1998) (see their Table 14), in 
which the lifetimes are calculated as the sum of the H-buming 
and He-buming timescales of the stellar tracks of Padua li- 
brary. In our models, for each mass M and metallicity Z 
we calculate the corresponding lifetime tmz by interpolating 
within the logarithmic relation log(M) - log(T) for the tabu- 
lated metallicities, and then by interpolating with respect to 
Z. 

2.5. Type la supernova 

Type la supernovas (SNe la) are important contributors of 
heavy elements, especially iron. Here we adopted the sin- 
gle degenerate model for the progenitors of SNe la, namely a 
degenerate C-O WD plus a red giant or main-sequence com- 
panion. The rate of SNe la is: 



Mb. 



RsNebit) - A 



Mb„, 



fifj) lp{t - Tmj) dju 



dniB (5) 



where A is a free parameter representing the fraction in the 
IMF of binary systems with the right properties to give rise 
to SNe la. Mbm and Mbm are the upper and lower limit re- 
spectively for the total mass of a binary system able to pro- 
duce a SNe la. As we use the chemical ejecta of stars by 



models with convective overshoot (see Section 4), the cor- 
responding upper limit to the mass of the primary is 6 Mq, 
so Mbm =12 Mq, while M^,,, ~ 3 Mq is adopted generally 
(Portinari et al. 1998). The function /(//) describes the dis- 
tribution of the mass ratio of the secondary (yu - M2/Mb) of 
the binary systems, and f(fi) - 2^^'''(l + y) jj^ with y = 2 as 
adopted by Greggio & Renzini (1983a). 

In this paper, we assum e the ejecta of SNe la are released 
in two steps as adopted in IPortinari et al.l (119981) : at the end 
of the primary's life, it expels its products just like a single 
star, while after the period of the lifetime of the companion 
the supernova explodes. So the contributions of SNe la are 
calculated in two parts: 
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where the first term in the square bracket represents the pri- 
mary's contribution and the second term represents the con- 
tribution of the SNe la explosion. 

2.6. Radial gas flow 

Radial gas flow in Galact ic disk can be driven by several 
mechanisms as introduced in|Lacey & Fall (1985), such as the 
difference between the angular momentum of the infalling gas 
and that of the circular motions in the disk, the viscosity in 
the gas layer, and the interactions between gas and the spiral 
structure. As the radial gas flows are expected to be inflows 
over most of the disk (Portinari & Chiosi 2000), we only con- 
sider inflows in our model, a nd we describe the g as infl ows 
following the prescriptions in lSpitoni & Mat teucci|(l201l1) . 

Let the k-th ring be represented by the galactocentric radius 
rii with its inner and outer edge being r^_i and r^^i (r^ i - 

(fk-\ + rk)l2 , = {rk + rk+i)/2). The velocities of the gas 
flowing through these edges are t;j_i and y^^i respectively 
with positive outward and negative inward. The surface gas 
density of the k-th ring is changed due to gas flows with a flux 
F{r) in the following way: 
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where the gas flow at r^^i can be written as: 
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Then we get: 
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where and jk are : 



a-1 + n 



n + n+i 
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We assume there are gas flows from the outer parts of the 
disk where is beyond r = 18 kpc, although there is no star 
formation in our model. We calculate the amount of gas flow- 
ing from the ring at r = 19 kpc to the ring at r = 18 kpc 
assuming that the surface gas density of the ring at r = 19 kpc 
evolves just like other rings which have star formation but its 
chemical composition stays primordial. This method may be 
not reasonable enough but it can easily help to reproduce the 
observational distribution of the gas along the radius at outer 
disk. Besides, the purpose of this paper is to study the effects 
of radial stellar migration, so we haven't spent much energy 
on the gas flow part. More detailed studies about the gas flow 
can be found in previous papers (e.g. iPortinari & Chiosil2000l : 
[Spitoni & Matteucci 201 1). 

The velocity of the gas flows adopted here is variable with 
radius as according to our tests, variable velocity is more help- 
ful to reconstruct the radial abundance gradients along the 
Galactic disk, and is also more reasonable. 



As the effects of all the factors above have already been 
studied in detail in many previous paper s (e.g. Portinari et al. 
[1991 iPortinari & Chiosil l2000t ISpitoni & Matteucci .2011.) . 
and the purpose of this paper is to find out the influence of 
the radial stellar migration on the Galactic chemical evolu- 
tion, we just adopt some typical prescriptions for the parts 
except stellar migration in the models as the bases to study 
stellar migration. The main parameters adopted for all the 
models are listed in Table [T] Two kinds of Galactic disk for- 
mation scenarios are adopted: 1. the whole disk formed si- 
multaneously with a constant disk formation time scale for 
all the Galactocentric radii, and these models are titled with 
'GCMl'; 2. the disk formed following an inside-out forma- 
tion scenario with the disk formation time scale increasing 
with radius, and these models are titled with 'GCM2'. The 
'(e)' at the end of the name of a model including radial stel- 
lar migration indicates that the values of the parameters on 
the right are adopted for both of the models adopting standard 
and extreme condition respectively for stellar migration (see 
Section 3), eg. GCMIN and GCMlNe. GCMl and GCM2 
are the basic models without radial stellar migration. Figure[T] 
shows the time-dependent radial iron and oxygen abundance 
distributions calculated by model GCM2 as an example. 

3. DESCRIPTION OF THE RADIAL STELLAR MIGRATION 

We describe the radial stellar migration in an analytical 
way, using different distribution functions to model the proba- 
bility that stars, bom at a given radius, appear at different radii 
after migrating for a period of time. This method doesn't refer 
to any single star but describes the radius distribution of stars 
quantitatively in a statistical way. Different from single stars, 
the migration of binaries may be more complex, but as the bi- 
naries always orbit around each other closely, we just assume 
the binary systems migrate like single stars. Thus the progen- 
itors of SNe la would migrate radially, the effects of which are 
considered in our calculation of the chemical contributions of 
SNe la. 
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Fig. 1 . — Galactic radial iron (top) and oxygen (bottom) abundance distribu- 
tion calculated by model GCM2, the lines from the bottom to the top repre- 
sent the abundance distributions from f = 4 Gijr to ( = 20 Gyr, respectively. 

The distribution functions we adopt for radial stellar migra- 
tion are normal distribution with stellar-age-dependent cr and 
Gama distribution with stellar-age-dependent shape and scale 
parameters. For the normal distribution. 



fsM = fN(r) = 
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where fi is the expectation corresponding to the maximum 
and the center of fNif), and cr is the variance representing the 
width or the shape of the distribution. We make the center at 
the birth radii of stars and cr increase with the stellar age expo- 
nentially so that stars would stay around their birth radii with 
the highest probability and migrate further with time elapsing 
(see the panel for birth radius R - 9 kpc in Figure |2]l. The 
normal distributions for all radii are adopted t o be same. 

Foll owing the hint of Figure 13 of Sell wood & Binneyl 
(12002b which shows the distributions of radii of stars with 
different initial birth radii, we adopt Gama distribution r{a,/3) 
which is a kind of asymmetrical distribution. The basic form 
ofr(Q',y8) is: 



fsM = Mr) = 
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exp -- , r > 



(10) 



r(a) 



CO 
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where a is the shape parameter and /3 is the scale parameter 
Considering the shape of Gama distribution will become more 
symmetrical as the shape parameter a increasing and wider 
with the scale parameter /3 getting bigger, we make a radius- 
dependent indicating stars bom at different radii would have 
different migration tendencies or preferences, and /3 would in- 
crease with the stellar age exponentially. At the same time, 
we also make the maximum of the Gama distributions corre- 
spond to the birth radii of stars. At 7? = 9 kpc, where is a 
transition radius, we adopt the normal distribution as same as 
in the pure normal distribution introduced above. 

For both the normal and Gama distributions, we adopt two 
scenarios for the variation of their profiles with stellar age: 
standard condition and extreme condition. In the extreme 
condition, stars would migrate much faster than in the stan- 
dard condition, which can't be possible in reahty. We just 
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TABLE 1 
Main parameters of the models. 



Models Stellar migration To ( A v "gfi^) 

(Gyr) (Gijf-^ ) {km s'^) 

GCMl No 3 0.25 0.07 0.55 -0.08r - 0.2, r < 10 ytpc 

GCMlN(e) Normal 3 0.25 0.07 0.55 -0.2r + 1.0, 10 /tpc < r < 15 <:pc 

GCMlG(e) Gama 3 0.25 0.07 0.55 -l/3r + 3.0, r > 15 ipc 



GCM2 No 0.88r-0.9 0.3 0.09 0.9 -0.16r + 0.6, r < 10 ytpc 

GCM2N(e) Normal 0.88r-0.9 0.3 0.09 0.9 -0.2r + 1.0, 10 /tpc < r < 15 /tpc 

GCM2G(e) Gama 0.88r-0.9 0.3 0.09 0.9 -\ /3r + 3.0, r > 15 kpc 

Note. — The '(e)' at the end of the name of a model including radial stellar migration indicates that the values of the parameters on the light are adopted for 
both of the models adopting standard and extreme condition respectively for stellar migration, eg. GCMIN and GCMlNe. To is the disk formation time scale, 
f is the fraction of the total stellar mass distributed in stars above 1 Mq, A is the parameter representing the fraction in the IMF of binary systems with the right 
properties to give rise to SNe la, v is the efficiency of the star formation rate, Vgf(r) is the velocity of radial gas flow. 



TABLE 2 

Values of ct in the normal distribution, a and /3 at all radii in the Gama 

DISTRIBUTION IN BOTH STANDARD AND EXTREME CONDITIONS. 



Normal 


Radius 








CT 


(kpc) 
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extreme 


allR 




exp(0.088f) - 


0.56 


exp(0.15t) - 0.56 


Gama 


Radius 


a 


/? 






(kpc) 




standard 




extreme 


4 
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£.1-/7(0.071/') - 


0.72 


exp(0.5t) - 0.12 


5 


4 


expiOMlf) - 


0.77 


exp{0.5t) - 0.77 


6 


5 


ejcp(0.054f) - 


0.80 


e,tp(0.5f) - 0.80 


7 


6 


exp(0.0495t) - 


-0.82 


exp(0.5f) - 0.82 


8 


10 


exp(0.038f) - 


0.86 


exp{0.5t) - 0.86 
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Normal 








10,11 


-10 


exp(0.038r) - 


0.86 


exp(0.5t) - 0.86 


12,13 


-6 


e.rp(0.04950 - 


- 0.82 


exp(0.5t) - 0.82 


14,15 


-5 


exp(0.054t) - 


0.80 


exp(O.St) - 0.80 


16,17,18 


-4 


exp{0.06U) - 


0.77 


exp{Q.5t) ~ 0.77 



Note. — t represents the age of star In Gama distribution, for the transition 
radius R = 9 kpc, the noimal distribution as same as in the pure normal 
distribution is adopted. The negative value for a represents we adopt the 
symmetrical distribution of the function with the positive value of the a. 

want to test what the effects would be if stars migrate under 
an extreme scenario. Table |2] shows the values of cr in nor- 
mal distribution, and a and at all radii in Gama distribution 
in both standard and extreme conditions, for the transition ra- 
dius R - 9 kpc the normal distribution as same as in the pure 
normal distribution is adopted. In Figure |2] we show the pro- 
files of Gama distribution at different radii for stellar age of 
2 Gyr (taller line) and 6 Gyr (lower line), and the normal dis- 
tribution at all radii is the same as that at 7? = 9 kpc. The solid 
line and dash-dotted line represent the standard and extreme 
condition, respectively. Here we point out that for simplicity 
the migration probability adopted here is constant in time. 

After migration, stars that die at radii different from their 
birth radii would eject their materials, including the newly 



synthesized materials and the ones taken from their birth radii, 
into the ISM at the positions they died, to some extent causing 
the mixture of material in large radius range in the disk. This 
may exert some effect on the Galactic chemical evolution. 

4. CONTRIBUTION OF STELLAR EJECTA AND 
NUCLEOSYNTHESIS PRESCRIPTIONS 

Stellar ejecta or yields are of significant importance in 
chemical evolution models. As the chemical enrichment of 
galaxies is basically due to stars, in spite of the light elements 
produced by primordial nucleosynthesis, the ejected materials 
of stars decide the amount of chemical elements in the ISM to 
a great extent. Here we adopt the stellar ejecta and yields cal- 
culated by models w ith convective overshoot (Marigo 20011 
iPortinari et al.lll998h for which the mass grouping is different 
from that of stand ard models (e.g. Woosley & Weaver 19861 
IChiosi et al.lll992h . with lower critical masses. Here we only 
briefly show the definition of different mass groups which are 
classified according to the dominant physical processes gov- 
erning their evolution and fate. For low-mass stars the upper 
mass limit is about 1.6 ~ 1.7 Mq in models with convec- 
tive overshoot, and about 6 Mq for intermediate-mass stars, 
both depending on the chemical composition; quasi-massive 
stars will explode as "electron capture" supernovas, the mass 
range of which is about >~ 6 to 8 Mp; and massive s t ars are 
in the mass range 8 ~ 120 Mq. (See IPortinari et alj (11998b 
for more details. Here we changed their upper mass limit of 
intermediate-mass stars from 5 Mq to 6 Mq, in consideration 
of the detailed explanation in their article.) 

For low- a nd intermediate -mass stars (m < 6 Mq) we adopt 
the yields bv lMarigol ( 1200 Ih which calculates the stellar evo- 
lution tracks up to the TP-AGB and III dredge-up phase, and 
the ejecta of massive and quasi-massive stars (6 ~ 120 M©) by 
IPortinari et al.l(ll998h . where the mass loss by stellar wind has 
been considered for massive stars, are adopted here. These 
yields and ejecta are metallicity dependent and calculated bas- 
ing on the same stellar models with the same series of stellar 
tracks for the whole range of stellar masses, so they are of 
good coherence and homogeneity of basic physical prescrip- 
tions. Fo r type la SNe w e adopt the ejecta EsNeia of the W7 
model of llwamoto et alJ ([1999). We assume each single star 
expels its ejecta all at once at the end of its life, while the bi- 
naries with the right properties to give rise to SNe la eject in 
two steps as introduced in Section 2.5, and the ejected mate- 
rial is immediately mixed in the ISM which remains always 



6 



0.8 \ 



0.6 ^ 



o 0.4 - 



0.2 



0.0 







2 4 6 



10 t2 2 4 



10 12 4 6 



10 12 124 4 6 



10 12 14 



0.8 



0.6 ^ 



o 0.4 - 



0.2 - 



0.0 



10 12 
r(kpc) 



1 4 m 



10 

r(kpc) 



12 



14 165 8 



10 12 14 
r(kpc) 



10 12 14 16 18 20 
r(kpc) 



Fig. 2. — Gama distribution for radial stellar migration for stars with birth radii /? = 4, 6, 8, 9 kpc (first row from left to right), and R = 10, 12, 14, 16 kpc (second 
row from left to right). The solid line and dash-dotted line represent the standard and extreme condition, respectively. For the same kind of line, taller lines are 
for 2 Gyr-old stars and lower lines are for 6 Gyr-old stars. The normal distribution at all radii in the pure normal distribution is the same as that at i? = 9 kpc. 



homogeneous. 



5. OBSERVATIONAL DATA 



of ir on we u se 
' ('2002c"a'.'b', 
(2011), 



For the Galactic abundanc e distribution 
the following sets of data: Andrievskv et al 
|2004), Lucketal. (2003, 2006), Luck & Lambert 
and ILemasle et al.l (120071, I2008D . who analyzed the Galac- 
tic Cepheids. For the abundance of oxyge n we adopt the 
data sets: lAndrievsk v et al.l (l2002ciyiR |2004^ . ILemasle et al.l 
(I2007h . and .LucketalJ (120031. 120061). who analyzed 
the Ci a lactic Cepheids; [S haver et al. ( 1983), Esteban et al. 
j2005'),'Vilchez & Estebml (11996), Rudolph et all (11997 ), and 
Rudolph et al. (2006), who studied the Galactic HII regions; 
Smartt & Rolleston ( 1997j), iGummersbach et al.1 (l998), and 
Daflon & Cunha (2004) , who analyzed the Galactic OB stars. 



Considering the scatter of the data, we divided the data into 
6 bins as functions of galactocentric distance, and calculated 
the mean abundance value and star-number-weighted radial 
position for each bin, to see the trend of the data more clearly. 
The mean value trends are represented by dark orange lines 
with filled circles in Figure [TOlandfTTI 

The data for the surface gas density profile along the Galac- 
tic disk are adopted from Dame (1993). Concerning the 
metallicity distribution of the G-dwarf stars in the solar neigh- 
borhood, we use the data of Rocha-Pinto & Maciel (1996)) . 
lJ0rgensenl (l2000h and lKotoneva et alj ( l2002h . 

6. RESULTS 

In this section we present the results of our models. 

In Figure |3] we show the surface mass density of gas pro- 



duced by our models. Solid lines and dash-dotted lines rep- 
resent the results of models with stellar migration of stan- 
dard and extreme condition, respectively. Dashed lines show 
the results of basic models without stellar migration. Lon; 
dashed lines show the observational data from iDamel (1 1 9931 
Comparing to basic models, there is less gas in regions where 
stars are lost by migration, since the materials ejected by 
these stars are also lost. As in the extreme condition stars 
migrate radially very fast, the gas produced by models with 
extreme stellar migration are less than that produced by mod- 
els with standard stellar migration, especially in the inner 
disk where more stars are born. However, the influence of 
the radial stellar migration on surface mass density of gas 
is small. Figure |4] shows the surface stellar mass density in 
the solar neighbourhood evolving with time. For the stan- 
dard condition, as more stars migrate out than into the so- 
lar vicinity, the stars produced by the models with normal- 
distribution stellar migration are less than that produced by 
basic models. But the stars produced by the models adopt- 
ing Gama-distribution stellar migration are more than that 
produced by the basic models. This results from that more 
stars born in the inner disk would migrate to the solar neigh- 
borhood under the Gama-distribution radial stellar migration. 
However, the current surface stellar mass density in the solar 
neighborhood produced by all the models with standard stel- 
lar migration can roughly fit the obser vational values, such 
as 63 Mn PQ-^ bvlJusi & Jahreifil (120101) and 61 ± 6 Mp pc' ^ 
bv lHolmberg & Flvnnl (l200'or and lHoknberg & Flvnnl dSool . 
Since in those papers only the part under vertical height 
z = 1 . 1 kpc above the disk plane was taken into consideration. 
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Fig. 3. — Surface mass density distribution of gas along the Galactic radius. 
Solid lines and dash-dotted lines in each panel represent the results of mod- 
els with stellar migration of standard and extreme condition, respectively. 
Dashed hues represent th e results of b asic models without stellar migration. 
Observational data from l5ani3 < 19931) are represented by long -dashed lines. 
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Fig. 4. — Surface stellar mass density in the solar neighbourhood evolv- 
ing with time. The solid lines and dash-dotted lines in each panel represent 
the results of models with stellar migration of standard and extreme condi- 
tion, respectively. The colors of the lines, blue and red, represent the results 
of models adopting normal and Gama distributions for the stellar migration, 
respectively. The black dashed lines represent the results of basic models 
without stellar migration. 

these values should be higher for the whole disk column. Thus 
the results calculated by the models with Gama-distribution 
radial stellar migration of standard condition should be closer 
to the real value. However, under the extreme stellar migra- 
tion, large amount of stars migrate out of the solar vicinity 
very quickly. Thus the surface stellar mass density produced 
by models with extreme stellar migration show a trend of first 
increasing and then decreasing, and the current stars are much 
less than the observational value. 

Considering the radial stellar migration, the G-dwarf metal- 
licity distribution of the solar neighborhood should be differ- 
ent from the results from star-static models (see Figure |5). 
As stars migrate radially, the solar neighborhood is not only 
assembled from local stars following a local age-metallicity 
relation, but also from stars originating from the inner (more 
metal-rich) and outer (more metal-poor) Galactic disk which 
have mi grated to the current position with different timescales 
dRoskar et al...2008: .Schonrich & Binnev.2009.) . Because of 



Fig. 5. — G-dwarf metallicity relation in the solar neighborhood by models 
with constant disk formation time scale (left) and inside-out disk formation 
scenario (right). In each figure, the dark orange line represents the result 
from the basic model without stellar migration; the blue and red lines rep- 
resent the results from the models including standard stellar migration with 
normal and Gama distribution respectively; the brown and green lines rep- 
resent the results from the models including extreme stellar migration with 
normal and Gama distribution respectively. The dashed, dash-dot ted and 
long-dashed histograms represent tlie observational data from Kotone va et alj 
f2OO2) , ,j0rgensen ( 2000) and .Rocha-Pinto & Maciel 1 ,1996) . respectively. 

the higher density of stars in the inner disk, the migrated 
portion would favor metal-rich stars. As a consequence, 
the migrated metal-rich stars could compensate the metal- 
poor tail of the local chemical evolution, and the profiles of 
the G-dwarf metallicity distributions we derived become nar- 
rower and taller with less metal-poor stars, while the mod- 
els with extreme condition really produce the extreme results, 
much more metal-rich stars. Thus, the stellar migration can 
be another solution to the "G-dwarf problem". However, 
the models adopting an inside-out disk formation, GCM2N 
and GCM2G, produce too many metal-rich stars in the so- 
lar neighbourhood, even though models adopting this kind of 
disk formation could have given out better results when do 
not consider the stellar migration. Thus radial stellar migra- 
tion may become another mechanism for the disk formation 
competing with the inside-out disk formation, or they can co- 
exist having similar effects. 

In Figure |6] we show the distributions of birth radii of the 
stars that appear at the solar neighbourhood at present. As 
we already know, the stars in the solar neighbourhood could 
have been born in a large range of radius and migrated to their 
present positions. While most stars come from the local area, 
a large part of the stars come from inner disk, as a conse- 
quence of the higher density of stars there, which is more 
obvious for the models adopting Gama distribution for stel- 
lar migration in which stars in inner disk have a preference to 
migrate outwards. These phenomenon are more obvious for 
the extreme condition. For our best model GCM2G, we show 
the distribution of the birth radii of stars in different ages that 
appear at the solar neighbourhood at present in Figure We 
can see that with the age increasing, the birth radii of stars 
cover larger range of the disk, and more stars come from in- 
ner disk, which is consistent with the stellar migration theory. 
We also show the time-dependent variation of the birth radii 
distributions for stars in the solar vicinity from the Galactic 
evolution time f = 1 Gyr to f = 14 Gyr in Figure |8] As the 
time goes, stars migrate for longer time, and more stars born 
at farther places come to the solar vicinity. Particularly for the 
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Fig. 6. — Distribution of the birth radii of stars that appear at the solar neigh- 
bourhood at present. Solid and dash-dotted lines in each panel represent the 
results of models with stellar migration of standard and extreme condition, 
respectively. The colors of the lines, blue and red, represent the results of 
models adopting normal and Gama distributions for stellar migration, respec- 
tively. 
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Fig. 7. — Distributions of the birth radii of stars in dilferent ages that appear 
at the solar neighbourhood at present by model GCM2G. Blue line is for stars 
younger than 3 Gyr; red line for stars of 3-7 Gyr; dark orange line for stars of 
7-11 Gyr; and lilac line for stars older than 1 1 Gyr. 

Stars from inner disk, we think that after about f = 4 Gyr they 
significantly enrich the stellar abundance in the solar vicin- 
ity, and this is also an explanation for the change of the G- 
dwarf metallicity relation caused by stellar migration which 
is shown in Figure |5] 

The age-metallicity relation (AMR) of stars in the solar 
neighbourhood is an important factor when considering stel- 
lar migration. As a consequence of radial stellar migration, 
the solar neighbourhood consists of stars of different ages 
and metallicities coming from different radii, and they scatter 
the classic AMR to an age-metallicity distribution (AMD), as 
roughly shown in Figure |9] The lines in each panel, from the 
top to the bottom, represent the contributions of stars migrat- 
ing from their birth radii of r = 4 kpc to 18 kpc respectively, 
and the relative proportion are shown by the colors (increasing 
from blue to red). As what we have expected, the AMR are 
scattered by the radial stellar migration, while a slope in the 
AMD is preserved and stars from the same birth radius still 
follow a classic local AMR. This phenomenon is also found 
bv .Minchev et al.. (.201Z) . 
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Fig. 8. — Time-dependent variation of the birth radii distributions for stars 
in the solar vicinity by model GCM2G, for the Galactic evolution time from 
t = 1 Gyr to f = 14 Gyr. 

Figure [TO] and [TT] show the comparison among abundance 
distributions from models and observation. The panels on the 
first and second row are for models with constant disk forma- 
tion time scale and inside-out disk formation scenario, respec- 
tively. From the left column, we can see that our models can 
reproduce the observational oxygen abundance and [Fe/H] 
distributions along the Galactic disk on the whole. In the 
right column we only compare the results from models. The 
abundance distributions in each panel, produced by the mod- 
els adopting the same Galactic disk formation scenario, dif- 
fer very little, almost coincide with each other Being incon- 
sistent with previous predictions, looking from the whole ra- 
dius range, the distributions of chemical element abundances 
along the disk have not flattened after radial stellar migra- 
tion, even in the extreme condition of radial stellar migration 
the abundance gradients have also not changed much. The 
specific values of the abundance gradient in the solar vicin- 
ity calculated by models are shown in Table |3] These gra- 
dients are a little flatter than the typical values which are 
obtained by fitting the whole data set along the disk, such 
as for the iron gradient: -0.070 by Lemasle et al. (2007); 
-0.052 by ,Lem asle et alJ (f2008); -0.051 by Pedicefli et all 
(120091) : -0.062 by Luck & Lambert (201 J), and for flie oxy- 
gen gradient: -0.044 by lEsteban et aL ('2 005h: - 0.06 by 
Rudolph et al. (2006) ; -0.06 5 by Lemasle et alT EoOTl) ; -0.055 
by iLuck & Lambert! (1201 I h. However, if only focus on the 
local gradient in the solar vicinity, we can find the gradient 
val ues for iron: -0.0 29 by |Andrievskv et al. (2002c); -0.03 
by iLuck et alj (f2003h; -0.0 44 bv ,Andrievskv et al.. (2004); - 
0.042 by iLuck & Lambert! (1201 Ih . which are consistent with 
our model results. 

Then we may infer that the radial stellar migration can not 
lead to significant flattening of the chemical abundance gra- 
dients. This is not diflicult to understand. As the impor- 
tant chemical contributors are massive stars, which have a 
short lifetime, they contribute mostly around their places of 
birth. Those migrated stars, born from far away having big 
difference in metallicity with the local stars, only contribute a 
small fraction to the local ISM and cannot change the chemi- 
cal abundance of the local ISM essentially; and the chemical 
contributions by stars from nearby differ not much enough 
to change the abundance gradients significantly. However, as 
the migration is a motion of stars, its influence on the stel- 
lar chemical abundance is more significant. It can flatten the 
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Fig. 9. — Age-metallicity distribution in tlie solar neighbourliood. In eacli panel from the top to the bottom, the lines represent the stars migrated to the solar 
neighbourhood from their birth radii of r = 4 kpc to 18 kpc, respectively, and the relative proportion are shown by the colors increasing from blue to red. 



radial gradient of the mean chemical abundance of stars. As 
in Figure [12] we show the radial distributions of the mean 
iron and oxygen abundance of stars in four age groups by 
model GCM2 and GCM2G. Since older stars have migrated 
for more time, the influences of radial stellar migration on 
them are more significant, making their chemical abundances 
higher and abundance gradients flatter in the solar nei ghbor- 
hood, which is consistent with the findings of Minch ev et alj 
(120121) . However, in the outer disk, as the disk there has 
formed much later and stars have a tendency to migrate in- 
ward under Gama-distribution migration, the chemical en- 
richment is slower and the mean stellar abundance gets lower 
for younger stars. Thus young stars have steeper abundance 
gradients in model GCM2G than in model GCM2. Besides, 
the most significant effect of the migrated stars on the chemi- 
cal abundance at a position is they can scatter the abundance 
of stars there from a relatively concentrated value to a range. 
As shown in Figure [T3] the shade shows the amount of star 
with a certain abundance value, which is calculated by model 
GCM2G as an example, and as a comparison the solid lines 
show the chemical abundance of the ISM. We can see that 
although the abundances of stars spread in a range, the abun- 
dance values which most stars possess at each radius still co- 
incide with that of the ISM. This, from another way, suggests 
that the radial stellar migration cannot change the abundance 
gradient of the ISM essentially. Moreover, as the gradients 
here are not very steep, the effect of scatter becomes more 
significant. 

We also predict that the abundance gradients of the Milky 
Way, such as the gradients of oxygen and iron, will not flatten 



TABLE 3 

Abundance gradients of iron and oxygen in the solar vicinity. 



Model 


d[Fe/H]/dR 


d[log(0/H)+12]/dR 




{dex kpc^^ ) 


{dex kpc^^ ) 


GCMl 


-0.036 


-0.039 


GCMIN 


-0.033 


-0.041 


GCMlNe 


-0.036 


-0.041 


GCMIG 


-0.033 


-0.042 


GCMlGe 


-0.036 


-0.041 


GCM2 


-0.041 


-0.033 


GCM2N 


-0.040 


-0.035 


GCM2Ne 


-0.042 


-0.033 


GCM2G 


-0.040 


-0.036 


GCM2Ge 


-0.042 


-0.033 



under the influence of radial stellar migration in the future. 
On the contrary, while the gradients have kept steepening for 
the past 13.7 Gyr, they will continue steepening in the future, 
but with slower speed, as shown in Figure [T4l 

7. SUMMARY AND CONCLUSIONS 

The aim of this paper is to study the influence of radial stel- 
lar migration on the chemical evolution of the Milky Way, 
especially for the Galactic abundance gradients. We did our 
research mainly based on two detailed Cialactic chemical evo- 
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Fig. 10. — Oxygen abundance along the disk obtained from models and observation. Panels on the first and second row are for models with constant disk 
formation time scale and inside-out disk formation scenario, respectively. Panels on the left column show the observational data and the results of models with 
standard radial stellar migration. The rhombuses, triangles and squares represent the observational oxygen abundance of Cepheids, HIl regions and OB stars, 
respectively. The dark orange line with filled circles shows the mean value of the observational data. Panels on the right column only show the results of models, 
adding the results of models with extreme radial stellar migration, so that we can compare the results of dilferent models clearly. In each panel, the black solid line 
represents the result of basic model without stellar migration; the blue solid line and the red dash-dotted line represent the results of models including standard 
stellar migration with normal and Gama distiibution, respectively; the golden long-dashed line and the lilac dashed line represent the results of models including 
extreme stellar migration with normal and Gama distribution, respectively. 



lution models: one adopts a constant disk formation time scale 
for all the Galactocentric radii; the other adopts an inside -out 
disk formation in which the disk formation time scale increase 
with radius. According to the basic point, we found the op- 
timal forms for the radial gas flow and values for the param- 
eters in the two models respectively to make them reproduce 
the observational data well. Then we added the radial stellar 
migration to these two models, where we used the distribution 
functions to model the radial stellar migration, and tested its 
effects in standard and extreme conditions, respectively. From 
the new models we derived the elemental abundance distribu- 
tions along the Galactic radius and other important informa- 
tion. 

Our main conclusions are summarized as follows: 

• When stars migrate radially on the Galactic disk, the so- 
lar neighborhood is not only assembled from local stars, 
but also from stars originating from the inner and outer 
Galactic disk. As there are more stars in the inner disk 
which are of higher metallicity, more metal-rich stars 
would migrate to the solar neighbourhood. As a con- 
sequence, the profiles of the G-dwarf metallicity distri- 
butions produced by our models considering the stellar 
migration are narrower and taller with smaller metal- 
poor tails than those produced by the models without 
stellar migration. Thus, the radial stellar migration can 



be another solution to the "G-dwarf problem". 

As model GCM2 adopting an inside-out disk forma- 
tion can produce taller and narrower G-dwarf metallic- 
ity distribution than model GCMl adopting a constant 
disk formation time scale, it seems that the radial stellar 
migration and the inside-out disk formation have some 
similar effects in building up the G-dwarf metallicity 
distribution. So whether the Galactic disk has formed 
inside-out or how it has formed in an inside-out way is 
still a question. 

As stars could migrate radially, the stars in the solar 
neighbourhood could have been bom in a large range of 
radius on the Galactic disk. Since older stars have mi- 
grated for more time, their birth radii cover larger range 
of the disk than younger stars, with larger fraction of 
them coming from inner disk. On the whole, under our 
assumption of the stellar migration, most stars in the 
solar neighbourhood come from the local area, while 
a large fraction comes from the inner disk, and fewer 
come from the outer disk. 

The classic AMR is scattered by the stellar migration 
to an age-metallicity distribution, while a slope in the 
distribution is preserved and stars from the same birth 
radius still follow a classic local AMR. 
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Fig. 11. — Iron abundance along the disk obtained from models and observation. Panels on the first and second row are for models with constant disk fomiation 
time scale and inside-out disk formation scenario, respectively. Panels on the left column show the observational data and the results of models with standard 
radial stellar migration. The dark orange line with filled circles shows the mean value of the observational data. Panels on the right column only show the results 
of models, adding the results of models with extreme radial stellar migration, so that we can compare the resu lts from different models clearly. 

The types of the lines are of the same meaning in figure [TOl 
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Fig. 12. — Distributions of the mean iron (top) and oxygen (bottom) abun- 
dance of stars in different age groups calculated by model GCM2 (dashed 
lines) and GCM2G (solid lines). Orange lines represent stars younger than 
3 Gyr; red lines represent stars of 3-7 Gyr; blue lines repres ent stars of 7- 
1 1 Gyr; green lines represent stars older than 1 1 Gyr 



Fig. 13. — Galactic iron (top) and oxygen (bottom) abundance distributions 
of stars by model GCM2G. The amount of star with a certain abundance 
value is shown by shade. As a comparison, the black solid lines show the 
abundance distributions of the ISM. 



We compared the results of models including radial 
stellar migration with those obtained from models with- 
out migration. Looking from the whole radius range, 
we find that the distributions of chemical element abun- 
dances along the disk haven't changed much after ra- 
dial stellar migration. Different from previous predic- 



tion, we infer that the radial stellar migration cannot 
change the chemical abundance gradients of the ISM 
essentially. However, it can flatten the radial gradients 
of the mean chemical abundance of stars. Since older 
stars have migrated for more time, they possess flatter 
abundance gradients than younger stars. 
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Fig. 14. — Galactic radial iron (top) and oxygen (bottom) abundance dis- 
tribution calculated by model GCM2G. the lines from the bottom to the top 
represent the abundance distributions from f = 4 Gyr to f = 20 Gyr, respec- 
tively. 



value to a range, while keep the dominant abundance 
value consistent with the abundance of the ISM. 

• According to the model results, we predict that the 
chemical abundance gradients of the Milky Way, such 
as the gradients of oxygen and iron, will continue steep- 
ening in the future, but with slower speed. 

Finally, we should point out that we just use two basic 
chemical evolution models as a foundation to study the ef- 
fects of radial stellar migration. We are not sure our models 
can produce the best results to fit all the observational data. 
Only according to our models in this paper, we also can not 
decide what kind of the combination of disk formation sce- 
nario, gas flow, stellar migration and the related parameters 
can make one model reproduce the observational data best. 
But we believe that our conclusions can, to a certain extent, 
tell the influences of radial stellar migration on the Galactic 
chemical evolution, and more further detailed research is still 
needed. 



• The most significant eff'ect of radial stellar migration on 
the chemical abundance at a position is it scatters the 
abundance of stars there from a relatively concentrated 
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